UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL 

EQUATIONS 

A. M. DAVIE 

Abstract. We consider the stochastic differential equation 

dx{t) = dW{t) + f{t, x{t))dt, x(0) = xq 

for t > 0, where x{t) e W is a standard d-dimensional Brownian motion, and / is 
a bounded Borel function from [0, ex.) X R'^ M'' to K'*. We show that, for ahiiost aU 
Brownian paths W{t), there is a unique x{t) satisfying this equation. 



1. Introduction 

In this paper we consider the stochastic differential equation 

dx{t) = dW{t) + fit, x{t))dt, x{Q) = xo 

for t > 0, where x{t) G M^, W is a standard (i-dimensional Brownian motion, and / is a 
bounded Borel function from [0, oo) x R'^ ^ R'^ to M'^. Without loss of generality we suppose 
Xo = and then we can write the equation as 

(1) x{t) = W{t)+ f f{s,x{s))ds, t>0 

Jo 

It follows from a theorem of Veretennikov [3] that ([T]) has a unique strong solution, i.e. there 
is a unique process x{t), adapted to the filtration of the Brownian motion, satisfying ([1]). 
Veretennikov in fact proved this for a more general equation. Here we consider a different 
question, posed by N. V. Krylov |2j: we choose a Brownian path W and ask whether ([1]) 
has a unique solution for that particular path. The main result of this paper is the following 
affirmative answer: 

Theorem 1.1. For almost every Brownian path W , there is a unique continuous x : [0, oo) — > 
M'^ satisfying (Qp. 

This theorem can also be regarded as a uniqueness theorem for a random ODE: writing 
x{t) = W{t) + u{t), the theorem states that for almost all choices of W, the differential 
equation ^ = f{t, W{t) + u{t)) with u{0) = has a unique solution. 

In Section HI we give an application of this theorem to convergence of numerical approxi- 
mations to ([I]). 

Idea of proof of theorem. The theorem is trivial when / is Lipschitz in x, and the idea of 
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the proof is essentially to find some substitute for a Lipscliitz condition. The proof splits into 
two parts, the first (section [2]) being the derivation of an estimate which acts as a substitute 
for the Lipschitz condition, and the second (section [3]) being the application of this estimate 
to prove the theorem. We start with a reduction to a slightly simpler problem. 

A reduction. It will be convenient to suppose \ f{t,x)\ < 1 everywhere, which we can by 
scaling. Then it will suffice to prove uniqueness of a solution on [0,1], as we can then repeat 
to get uniqueness on [1,2] and so on. 

So we work on [0,1], let X be the space of continuous functions x : [0,1] Mf^ with 
x(0) = 0, and let Pw be the law of M*^- valued Brownian motion on [0,1], which can be 
regarded as a probability measure on X. Now we apply the Girsanov theorem (see [3]): 
define (j){x) = exp{ f{t, x{t))dx{t) — | /(t, x{t)Ydt}, which is well-defined for Pw almost 
all X G X, and define a measure on X by d/x = (pdPw- Then if a; G X is chosen at random 
with law /i, the path W E X defined by 

(2) W{t) = x{t) - [ f{s,x{s))ds 

Jo 

is a Brownian motion, i.e. W has law Pw- 

For a particular choice of x, and with W defined by ([2]), x will be the unique solution of 
([1]) provided the only solution of 

(3) u{t)= [ {f{s,x{s)+u{s))-f{s,x{s))}ds 

Jo 

in X is M = 0. So, to prove the theorem it suffices to show that, for ([3]) has no 

non-trivial solution, since for such x, with W defined by ([2]) no other x can satisfy ([2]). 

But /i is absolutely continuous w.r.t. Pw, so it suffices to show that, for Pyjz-a.a. x, has 
no non-trivial solution. In other words, it suffices to show that, if is a Brownian motion 
then with probability 1 there is no non-trivial solution m G X of 

(4) u{t) = [\f{s, W{s) + u{s)) - f{s, W{s))}ds 

Jo 

We prove this in section [3l 

Remark. Our proof does not make use of the existence of a strong solution. It is tempting 
to try to prove the theorem by measure-theoretic arguments based on the strong solution 
and Girsanov's theorem. Define T : X ^ X by 

Tx{t) = x{t) - [ f{s,x{s))ds 
Jo 

The strong solution gives a measurable map S : E —* F where E and F are Borel subsets 
of X with Pw{E) = Pw{F) = 1, such that T o is the identity on E, and F is the range 
of S. It follows that T is (1-1) on F and for any W E E there is a unique solution of ([1]) in 
F. But we need a solution which is unique in X and to achieve this we need to show that 
T(X\F) is a Pvi/-null set, and this seems to be a significant obstacle. 
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Our proof is quite complicated and it seems reasonable to hope that it can be simplified. In 
particularly one might expect a simpler proof of Proposition [521 This seems to be nontrivial 
even for p = 2. The bound for p = 2 follows from the first part of Lemma 12.51 (with to = 
and r = 0) and I do not know an essentially simpler proof. 

In one dimension, in the case when f{t, x) depends only on x, a different and shorter proof 
of Theorem II. II can be given, using local time, but it is not clear how to extend it to d > 1. 

2. The basic estimate 

This section is devoted to the proof of the following: 

Proposition 2.1. Let g be a B or el function on [0, 1] x M'^ with \g{s,z)\ < 1 everywhere. For 
any even positive integer p and x G M.'^, we have 

E ( [\git, Wit) +x)- git, Wit))}dt] < CPip/2)\\xf 



where C is an absolute constant, \x\ denotes the usual Euclidean norm and Wit) is a standard 
d- dimensional Brownian motion with WiQ) = 0, 

This will be deduced from the following one- dimensional version: 

Proposition 2.2. Let g be a compactly supported smooth function on [0,l]xM. with \gis, z)\ < 
1 everywhere and g' bounded (where the prime denotes differentiation w.r.t. the second vari- 
able). For any even positive integer p, we have 

-1 \P 



e(^^ g'it,Wit))dt^ <C\p/2)\ 



where C is an absolute constant, and here Wit) is one- dimensional Brownian motion with 

iy(o) = 0. 

Proof. We start by observing that the LHS can be written as 

p 

p 



E\[g'itj,Witj))dti---dtp 

0<ii<--<tp<l j^-^^ 



and using the joint distribution of W^(ti), ■ ■ ■ , Witp) this can be expressed as 

p\ I I Wig'itj, Zj)Eitj - tj_i, Zj - Zj-i)]dzi ■ ■ ■ dzpdti ■■■dtp 

Jo<ti<- -<tp<l Jrp j^-^ 

2 / 



where z) = (27rt)-i/2g-^ /2t ^lere to = 0, = 0. 
We introduce the notation 

r r ^ 



/ / W{9'itj,Zj)Eitj-tj-i,Zj- Zj^i)]dzi^^^dzkdti^^^dtk 
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and we shall show that Jp(0, 0) < C^/r(| + 1); Proposition 12.21 will then follow since 
p\ < 2f((p/2)!)2. 

In order to estimate Jk we use integration by parts to shift the derivatives to the ex- 
ponential terms. We introduce some notation to handle the resulting terms - we define 
B{t,z) = E'{t,z) and D{t,z) = E"{t,z) (where again primes denote differentiation w.r.t. 
the second variable). 

If S = Si ■ ■ ■ Sk is a. word in the alphabet {E, B, D} then we define 

k 

Is{to,zo)= / / Y\{g{tj,Zj)Sj{tj -tj_i,Zj - Zj^i)}dzi - ■ ■ dzkdti- ■ -dtk 

J tQ<tl<-<tk<l JVid 

In fact, only certain words in {E, B, D} will be required: we say a word is allowed if, when 
all -B's are removed from the word, a word of the form {EDY = EDED ■ ■ ■ ED, r > 0, 
is left. The allowed words of length k correspond to the subsets of {1, 2, • • • , k} having an 
even number of members (namely the set of positions occupied by E and D in the word). 
Hence the number of allowed words of length k is the number of such subsets of {1, 2, ■ ■ ■ ,k}, 
namely 2^~^. 

We shall show that 

2fc-l 

(5) Jkito, zo) = ^ ±lsu) (to, zo) 

i=i 

where each 5**^-'^ is an allowed word of length k (in fact each allowed word of length k appears 
exactly once in this sum, but we do not need this fact). The proof will then be completed 
by obtaining a bound for Is- 

We prove ([5]) by induction on k. So, assuming ([5]) for Jk, we have 



dti J g'{ti,Zi)E{ti -to,Zi - Zo)Jk{ti,Zi)dzi 
- dti g{ti,zi)B{ti -to,zi - zo)Jk{ti,zi)dzi 

Jtn J 

g(ti, zi)E(ti - to, zi - ZQ)Jk{ti, zi)dzi 



to 



Now we observe that, if S is an allowed string then I'g = —Ig where S is defined as BS* if 
S = ES* and as DS* if 5* = BS* (note that S is not an allowed string). Applying this to 
we find J'^ito, Zq) = Yl^j=i -F -^53(^0; ^0) aiid then we obtain 

■Jk+lito, Zo) = T -^BS^ (^0, ^0) ± Ieso{'^o-,zo) 



< c 
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Noting that, if S is an allowed string, BS and ES are also allowed, this completes the 
inductive proof of ([5]). 

We now proceed to the estimation of /^(to, zq), when S is an allowed string. We start with 
some preliminary lemmas. 

Lemma 2.3. There is a constant C such that, if (p and h are real-valued Borel functions on 
[0, 1] X M with \(l){t,y)\ < e'^ and \h(t^y)\ < 1 everywhere, then 

I dt ds / (l){s, z)h{t,y)D{t — s,y — z)dydz 

J 1/2 Jt/2 JrJr 

Proof. Denote the above integral by /. For / G Z, let xi be the characteristic function of 
the interval [1,1 + 1) and define (f)i{s,y) = (p{s,y)xi{y), and similarly hi. Let Iim denote the 
integral / with (p, h replaced by (pi, km- Then we have / = ^im- Let Ci, C2, ■ ■ ■ denote 

positive absolute constants. 

Now if |/ — m| = > 2 then for z E [/,/ + !) and y G [m, m + 1) we have \z — y\ > k — 1 
and then it follows easily that 

\D{t-s,y-z)\<Cie-^'-^^'^^ 
and hence Iim < (72e~'^/^e^''^^^^^/'^ from which we deduce 

\Ilm\ < C*3 

|«-m|>2 

Now suppose |/ — m| < 1. We use 4^i{s, u) for the Fourier transform in the second variable, 
and similarly hm- We note that / (pi{s,uYdu = J (pi{s,zydz < C4e~l'l for < s < 1 and 
similarly f hm(t,uYdu < 1. We have 

Iim= [ dt f ds [ ^i{s,u)h„^{t,-u)e''-'"^^''\'/\^du 

Jl/2 Jt/2 Jr 

Applying ab < |(a^c + b^c~^ with a = (pi{s,u), b = hm(t, —u) and c = e'^/"*^^, we deduce that 

\Iin.\ < f dt f ds [ Us,ufe'"/^We~^'-''>"'^du 
Ji/2 Jt/2 Jr 

+ f dt f ds [ h^{-t,ufe-'"/^We~^'-'>"'^du 
J 1/2 Jt/2 Jr 

In the first integral we integrate first w.r.t. t and obtain the bound const. e~'^/^^ for the 
integral. We get a similar bound for the second integral (integrating w.r.t. s first), and 
hence 

\hm\ < C,e-''l'' 
Summing over / and m such that |/ — m| < 1, we obtain 

\ilm\ ^ C*6 

lZ-m|<l 
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which completes the proof. □ 

Corollary 2.4. There is an absolute constant C such that if g and h are Borel functions on 
[0, 1] X R bounded by 1 everywhere then 



dt 

1/2 Jt/2 



ds I g{s, z)E{s, z)h{t,y)D{t — s,y — z)dydz 



and 



/ dt ds g{s, z)B{s, z)h{t,y)D{t — s,y — z)dydz 

Jl/2 Jt/2 Jr^ 



< c 



< c 



Proof. These follow easily from Lemma fl2.3p . the second using the easily verified fact that 

We note that f^E(t, z)dz = 1, and we have the bounds 

(6) [ \B{t,z)\dz <Cot-^/^, I \D{t, z)\dz < Cot-^ 

Jr Jr 

where Cq is an absolute constant. 

Lemma 2.5. There is an absolute constant C such that if g and h are Borel functions on 
[0, 1] X M bounded by 1 everywhere, and r > then 



to 

and 
1 



to 



dt ds g{s, z)E{s — to, z)h{t,y)D{t — s,y — z){l — tydydz 



to 



to 



dt ds g{s, z)B{s — to, z)h(t,y)D(t — s,y — z)(l — tydydz 



<C{l + r)-\l-toY+' 



< C(l+r)-i/^(l-to)''+^ 



Proof. Again, we let Ci, ■ ■ ■ be absolute constants. By using the change of variables t' = 
— to) /{I — to), s' = {s — to) /{I — to), y' = — to)~^^'^, it suffices to prove these estimates 
when to = 0. To do this, we start by scaling the first part of Corollary 12.41 and get 



dt 



2-k-l 



t/2 



ds / g{s,z)E{s,z)h{t,y)D{t — s,y — z){l — tydydz 



for = 0, 1, 2 ■ ■ ■ and then by summing over k, we get 



dt 



t/2 



ds / g{s, z)A{s, z)h{t,y)D{t — s,y — z){l — ty dydz 



Moreover, from the bounds ([6]) we have 



dt 



t/2 



< 



ds / g{s, z)E{s, z)h{t,y)D{t — s,y — z){l — ty dydz 
Jr^ 

j-l j-t/2 

Cs dt {t-s)-\l-tyds<Ci{l + r)-^ 
Jo Jo 



< 
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and combining these bounds gives the first result. Similarly, by scaling the second part of 
Corollary 12. 4^ we get 



dt 



2-fc-i 



t/2 



ds / g{s^ z)B{s, z)h{t,y)D{t — s,y — z){l — tydydz 



for /c = 0, 1, 2 ■ ■ ■ and then by summing over k, we get 



dt ds g{s^ z)B{s, z)h{t,y)D{t — s,y — z){l — tydydz 

Jo Jt/2 



't/ 

Moreover, from the bounds (El) we have 



<C.(l + r) 



-1/2 



/•I i-m r 

dt ds g{s, z)B{s, z)h{t,y)D{t — s,y — z){l — ty dydz 
Jo Jo Jm? 

<Co dt {t- s)-\l - tyds < Cril + r)-i/2 
Jo Jo 



< 



which give the second result. 



□ 



We can now complete the proof of Proposition 12.21 by obtaining the required bound for 
-^5(^05 zq). Again we use Ci, C2, ■ ■ ■ for absolute constants. We shall show that, for a suitable 
choice of M, we have for any allowed string S of length k 



(7) 



\Is{to,zo)\ < 



r(f + 1) 



{l-tof' 



We shall prove ([7]) by induction on k, provided M is chosen large enough. The case = is 
immediate, so assume k > and that ([7]) holds for all allowed strings of length less than k. 
Then there are three cases: (1) S = BS' where 5" has length — 1; (2) 5* = EDS' where 5" 
has length k — 2; {?>) S = EB"^DS' where m > 1 and S' has length k — m — 2. In each case 
S' is an allowed string. We consider the three cases separately. 

Case 1. In this case we have 

•1 



1/5(^0,2:0)1 



< 



dti / B{ti-to,zi- zo)g{ti,zi)Is'{ti,zi)dzi 

to Jk 

/ \B{t^-to,z,-zo)\dz, 



r 



■fe+1' 



to 



< 



CiM 



k-l /-l 



\l-t^f-'y\t^-to)-'/'dt^ 



to 



where we have used the inductive hypothesis to bound Is', and then the bound ([6]). ([7]) then 
follows if M is large enough. 
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Case 2. Now we have 

Is{to,Zo)= dti dt2 g{ti,zi)g{t2,Z2)E{ti-to,zi-zo)D{t2-ti,Z2-Zi)Is'{t2,Z2)dzidz2 

Jta Jti JR2 

We set h{t,z) = g{t, z)Is'it, z){l - ty-^ so that \\h\\oo < M''-^/T{k/2) by the inductive 
hypothesis, and then from the first part of Lemma 12.51 we deduce that 

\i(f ,^|<^i^^!:!(i^M!^ 
i/5(to,^o)i< — — 

and ([7j) follows if M is large enough. 
Case 3. In this case have 



„ /m+2 \ 

dti--- dtm+2 / TT 9{tj, Zj) E{ti - to, zi - zq) : 



Is{to, Zo) = dti--- dtm+2 I III 9{tj, Zj) I E{ti - to, Zi - Zo) X 

Jto<tl<---<tm+2<l 

m+1 

X W B(tj — tj^i, Zj — Zj^i)D(tm+2 ~ tm+l, Zm+2 ~ Z^+l) I S' {tm+2, Zm+2)dZi - - - dZm+2 
i=2 

Now let h{t,z) = g{t,z)Is'{t,z){l — so that by the inductive hypothesis on 

S' we have \\h\\^ < M^-'^-^/T{^). Then, writing 

Q,(t,z) = I dtm+l / dtM+2 / git-m+l, Zm+l)h{tm+2, Zm+2){^ — tm+2)^^ ™ ^''''^ X 

Jt Jtra+l JR'^ 

X B{tm+1 — t, Zm+1 — z)D(tfn+2 ~ tm+l, Zm+2 ~ Zm+l)dZm+ldZm+2 

we find from Lemma 2 that 

m, z)\ < Cs{k - m)-^'^M^-^-\l - t)('=-™-i)/2/r ~ 

Using this in 

h{to, Zo) -- 



dti--- dtm / n Zj) E{ti - to, zi - Zo) x 



'to<tl<---<tm<l 

m 

X Y[B{tj-tj _1, Zj Zj-i)Q{tm, Zm)dZi - - - dz,, 
i=2 

and using the bounds (jH]) we find 

j^k—m—2 

; Jto<ti<-<tm<l 

{tm tm—l) ^ (1 tm)^ ' dti - - - dtm 

]\/rk-m—2^(in-l)/2rf k—m+l \ 

- m)-'/'- - to)'/' 



Is{to,zo)\<Cr\k-m)-'/'^,yj-^ I 
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from which again ([71) follows, provided M is large enough. Putting ([7]) with to = 0, = 
and k = p in ^ completes the proof of Proposition 12.21 □ 

Proof of Proposition \2.1\ We first note that it suffices to prove it for d = 1. To see this 
let W, X be as in the statement of Proposition 12. 1[ By a rotation of coordinates we can 
suppose a; = (a, 0, ■ ■ ■ , 0). Then for fixed Brownian paths W2, ■ ■ ■ , Wd we can define h on 
[0, 1] X M by h(t, u) = g(t, u, W2(t), ■ ■ ■ , Wd(t)) and the d = 1 case of the Proposition gives 

e( [ {h{t,Wi{t) + a) - h{t,Wi{t))}dt] <CP{p/2)\\a\P 



and then the required result follows by averaging over W2, ■ ■ ■ , W^- 

So we suppose d = 1. Given a Borel function g on [0, 1] x M with \g\ < 1 we can find 
a sequence of compactly supprted smooth functions gn with \gn\ < 1, converging to g a.e. 
on [0,1] X R. Then gn(t,W(t)) g(t,W(t)) a.s. for a.a. t e [0,1], and the same for 
gn(t, W(t) + x), so by Fatou's lemma it suffices to prove the proposition for smooth g. But 
then we have g(t, W(t)+x) —g(t, W(t)) = g'{t, W(t) +u)du and we can apply Proposition 
12.21 and Minkowsi's inequality to conclude the proof of Proposition 12.11 

What we in fact need is a scaled version of Proposition 12.11 for subintervals of [0,1]. For 
s > we denote by J-'s the a-field generated by {W{t) : < r < s}. Then we can state the 
required result: 

Corollary 2.6. Let g be a Borel function on [0,1] x R"' with \g\ < 1 everywhere. Let 
0<s<a<b<land let p{x) = JI^{g{W{t) + x) - g{W{t))}dt. Then for x eR'^ andX>0 
wc hdvc 

nip(a;)l > A/'/'kl I < 2e-"'/(2c2) 
where I = b — a and C is the constant in Proposition \2.1[ 

Proof. First assume s = a = 0, b = 1. Let a = {2C'^\x\'^)^^ . Then 

00 ^ 00 
E(eM-)^) = J2^E{p{xf'') < ^a'^C'^lxp'^ = 2 

fc=0 ■ k=0 

and so P(|p(x)| > A|a;|) = P(e°''(^)' > e^^'l^'l') < 2e~"^'l^l' = 2e~^'/^'^^'\ 

For the general case, let W{t) = l^^^'^ {W (a+tl) —W (a)} , so that is a standard Brownian 
motion, and let h{x) = g{W{a) + x). Then p{x) = /^/^ J^{h{W{t) + x) - h{W{t))}dt and 
it follows from the first part that P(|p(x)| > A/^/^|x| < 2e^^^^^'^'-'^\ The required result 
then follows by taking conditional expectation w.r.t. JF,. □ 

We note that the unconditional bound 

F{\p{x)\>Xl'/^\x\)<2e-^"/^'^"^ 
follows by taking s = 0. Also in the same way we obtain, for any even p G M, 
(8) E(p(x)P|J-,) < CPF/\p/2y.\x\P 

The following lemma will also be needed: 
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Lemma 2.7. If p > 1 + f there is a constant c{p, d) such that is g & LP{\fi,l] x W^) then 

e( r g{t,W{t))dt) <c{p,d)\\gf 







p 



Proof. We have 

E f /" g{t,W{t))dt] =2 1 dt I ds I g{s,C)g{t, z)E{s,C)E{t - s, z - C)dCdz 



Now, if g = ^ then / E{t, zfdz = 0{t-^''-^^'^/^) and p > 1 + f imphes {q - l)d/2 < 1, so 
the result follows from Holder's inequality. □ 

3. Proof of Theorem 

We now apply Corollary 12.61 and Lemma 12.71 to the proof of the theorem. First we give a 
brief sketch of the proof. 

Outline of proof. The proof is motivated by the elementary case when / is Lipschitz in 
the second variable. In this case, if / = [a,b] is a subinterval of [0,1] and m is a solution of 
@ satisfying 

(9) \u{t)\<a, tel 

and P = \u{a)\, then we deduce from that \u{t)\ < a' = P + L\I\a for t E I, where L is 
the Lipschitz constant, i.e. holds with a replaced by a'. If L\I\ < 1 it follows that 
holds with a = (1 — L\I\)^^P, and of course if /? = this gives m = on /. 

We try to copy this argument using Corollary 12.61 as a substitute for a Lipschitz condition. 
There are two difficulties: first. Corollary 12. 6l is a statement about probabilities and we need 
an 'almost sure' version, and in doing so we lose something; second, in Corollary 12. 6[ x is 
a constant, whereas we are dealing with a function u depending on t. The way round the 
second problem is to approximate m by a sequence of step functions ui and then use 

(10) ^ 

'{fiWit) + u{t)) - fiWmdt = hm / {f{W{t) + ui{t)) - fiWmdt 



CO 



I l=n 



{f{W{t) + u^{t)) - f{Wmdt + 2^ / {f{W{t) + ui+m - f{W{t) + uimdt 



where u„ is constant on the interval /, and then to apply the 'almost sure' form of the 
proposition to each interval of constancy of the terms on the right. Again, we lose something 
in doing this, but, as it turns out, we still have good enough estimates to prove the theorem. 
In fact, we need two versions of the 'almost sure' (nearly) Lipschitz condition, the first to 
estimate / {f {W {t) + Un{t)) - f{W{t))}dt and the second to estimate / {f {W {t) + ui+i{t)) - 
f{W{t) + ui{t))}dt. We also need a third estimate, for sums of integrals of the second type. 

The two versions of the 'almost sure' nearly-Lipschitz condition are conditions f|TT]) and 
([12]) below, and the third estimate is (1201) . In Lemmas 13.11 13. 2[ 13.51 and 13.61 it is shown 
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that these conditions indeed hold almost surely. Lemmas 13.31 and 13.41 establish a technical 
condition f[T^ needed to justify the passage to the limit as / oo (which is not trivial 
when / is not continuous). With these preliminaries the above programme is carried out in 
Lemma [3.71 The analogue of (Q above is (l25ll . We no longer immediately get a = when 
P = 0, but we get a good enough bound to prove the uniqueness of the solution to ([1]), for 
any W satisfying (fTT][T2][T5]l20D . 

We now turn to the details. 

For any ri > we can divide [0,1] into 2" intervals Ink = [k2^'^'',{k + 1)2~"], k = 
0, 1, 2, ■ ■ ■ , 2" — 1. We shall also consider dyadic decompositions of M.'^, and say x G R'^ 
is a dyadic point if each component of x is rational with denominator a power of 2. Let 
Q = {x & M."^ : \\x\\ < 1}, where denotes the supremum norm maxi<j<d\xj\. We also 
introduce the notation 

c^nkix) = [ {g{W{t) +x)- g{W{t))}dt 

J Ink 

and 

Pnk{x, y) = (Tnkix) - ank{y) = / {g{W{t) +x) - g{W{t) + y)]dt 

''Ink 

Then we can state: 

Lemma 3.1. Let g he a real function on [0, 1] x with \g{t,z)\ < 1 everywhere. Then with 
probability 1 we can find C > so that 

(11) |p.fc(x,i/)| <c|ni/2+ (^log+^l^y^'|2-"/^|x-i/| 

for all dyadic x,y E Q and all choices of integers n, k with n > and < A; < 2" — 1. 

Proof. Let us say that two dyadic points x,y are dyadic neighbours if for some integer 
m > we have — y\\ = 2~™ and 2-™x, 2-"y G Z'^. Then using the Corollary 12.61 we have, 
for any such pair x,y E Q and any n, k that 

P {\pnk{x,y)\ > A(n^/2 ^m^/2)2— < (j^^-C,X^{n+m) 

and by summing over all possible choices of n, k, m, x, y we find that the probability that 

|p„fc(a;,2/)| >A(nV2+ ^1/2)2— "/2 

for some choice of Ink and dyadic neighbours x,y E Q is not more than 
Er=i ^;;^^o2"3'^2'^('"+3)Cie-^2A2(i+m+n) ^j^jpj^ approaches as A ^ oo. 

It follows that, given e > 0, we can find A(e) such that, with probability > 1 — e, we have 

|p„fc(a;,y)| < A(l + ni/2 + mi/2)2— n/2 

for all choices of n, k and dyadic neighbours in Q. 

Next let X, y be any two dyadic points in Q, and let m be the smallest non-negative integer 
such that \\x — y\\ < 2^™. For r > m, choose Xr to minimise \\x — Xr\\ subject to 2^Xr G Z'^, 
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and i/r similarly. Then \\xm — ym\\ = or 0, and for r > m, \\xr — Xr+i\\ = 2^^^-^ or 0. So 
Xm,ym are dyadic neighbours or equal, and the same applies to Xr,Xr+i and Ur^Ur+i- Then 
we have 

oo oo 

r=m T=m 

(note that the sums are actually finite, since x, y are dyadic, so that x = Xr and y = yr for 
large r). Then applying the above bounds for the case of dyadic neighbours to each term, 
we get the desired result. □ 

Next we prove a similar estimate for ank, which is analogous to the Law of the Iterated 
Logarithm for Brownian motion. 

Lemma 3.2. With probability 1 there is a constant C > such that for all n E 'H, k & 
{0, 1, ■ • ■ , 2" — 1} and dyadic x E Q we have 

(12) \a^k{x)\<Cn'/'2-^^/\\x\+2~'") 

Proof. For any integer r > we let = {x E M.'^ : \\x\\ < 2"*"}. Then if m > r the number 
of pairs (x, y) of dyadic neighbours in Qr with ||a; — y\\ = 2"™" is < (9 x 2™-"^)'^ and for each 
such pair we have 

F{\pnk{x,y)\ > A(nl/2 + V^^r^)2-'"~"/2) < ^^^^C2XHn+m-r) ^ (j^^2d{r~ra) ^~C-,\\-n 

for A large. By summing over n, 1 < r < 2" and m > r and all pairs {x, y), we deduce that, 
with probability > 1 — C^e''-''^'^ , we have Pnk{x,y) < Xip}/"^ + v^rn^^)2~^~"/^ for G N, 
1 < r < and m > r and all pairs {x,y) of dyadic neighbours in Qr with ||x — y\\ = 2"™", 
and then, by an argument similar to Lemma 13. H we get for all n and 1 < r < n that 
o'nkix) < C5An^/^2~^~"/^ for all dyadic x E Qr- The required result follows. □ 

The next two lemmas are used to justify the passage to the limit / — oo in (ITO!) . 

Let $ denote the set of Q- valued functions u on [0,1] satisfying \u{s) — u{t)\ < \s — t\, 
s,t E [0, 1], and let denote the set of Q- valued functions on [0,1] which are constant on 
each 4fc and satisfy |M(fc2-") - m(/2-")| < \k - /|2-". Then let = $ U U„$„. 

Lemma 3.3. Given e > 0, we can find r] > such that if U G (0, 1) x is open with 
\U\ < rj, then, with probability > 1 — e, we have Xuitj W{t) + u{t))dt < e for all u G $*. 

Proof. Fix e > 0. By Lemma 13.11 we can find K such that, for any Borel function (p on 
[0, 1] X M'^ with 101 < 1 everywhere we have with probability > 1 — e/2 that 

(13) f {(l){W{t)+x)-(l){W{t)+y)}dt<Kn'/^2-'^''/^ 

J Ikn 

for all pairs of dyadic points x,y in Q and all choices of n, k. Then we choose m such that 
AK Yl'^=m ^^^''^2~"'/^ < e. Let f2 be a finite set of dyadic points of Q such that every x E Q 
is within distance 2~™ of some point of VL. 
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Provided 6 is chosen small enough, any bounded Borel function on [0, 1] x M"* with 
LP([o,i]xR'i) < ^ will satisfy 



P 



W{t) + x)dt 



— ' I nm+l 



for each fc, x. Then the probability that 



(14) 



0(t, W{t) + x)dt 



< 2 "^6/4 for every k G {0, 1, 



1}, X E^l 



is at least 1 — e/2. 

Now let 7] = 6^, and suppose U is open with m{U) < rj. Let (0^) be an increasing sequence 
of continuous non-negative functions on [0, 1] x W^, converging pointwise to xu- Note that 
then ||0r||LP([o,i]xM<* < ^- For sach r define events A^: f|T^ holds for (p = (p^ and 5,.: f fT3l) 
holds for (p = '(pr. Then F{Ar) > 1 - e/2 and P(5^) > 1 - e/2. Also, when Ar and 5^ both 
hold, we have j]^^ (pr{t, W{t) + x)dt < 2"™e/2 for all x such that |a;| < 2. 

Now let e $*. For each n > m choose u„ G $„ taking a constant dyadic value within 
2-" of m(A;2-") on 4^ for A; = 0, 1, ■ ■ ■ , 2" - 1. Now if Ar and 5^ hold then (pr{t, W{t) + 
Um{t))dt < e/2 and 



{0,,(t, W{t) + u„(t)) - 0,(t, W{t) + M„+i(t))}dt 



< irn^/22-"/2 



from which it follows that (pr{t,W{i) + u{t))dt < e. So if we define the event Qr ■ 

Jq (prit, W{t) +u{t))dt < e for all u G (p, then we have F{Qr) > 1 — e. But since 0r+i > (pr we 
have Qr+i ^ Qr, and it follows that with probability > 1 — e we have Qr for all r, from which 
the result follows, since /q (prit, W{t) + u{t))dt /q + u{t))dt by the bounded 

convergence theorem. □ 

Lemma 3.4. Ifg is a bounded Borel function on [0, IjxM'^, then, with probability 1, whenever 
(m„) is a sequence in $* converging pointwise to a limit -u G $*, we have 

(15) [ g{t,W{t)+Un{t))dt^ [ g{t,W{t)+u{t))dt 

Jo Jo 

Proof. Given e > 0, let be as in Lemma [X^ and let hhe a. bounded continuous function on 
[0, 1] xR'^ such that g = h outside an open set U with m{U) < rj. With probability > 1— e, the 
conclusion of Lemma [3.31 holds, which means that for any convergent sequence (u„) in $ we 
have luit, W{t) +Un{t))dt < e, and the same for the limit u{t), so, if M is an upper bound 
for \g-h\, we have the bound \J{g{t,W{t) +M„(t)) - h{t,W{t) +M„(t))}dt| < Me, and the 
same for x in place of u„. Also, since h is continuous, J^ h{t, W{t) +Un{t))dt J^ h{t, W{t) + 
u{t)dt. It follows that, for n large enough, J^ g{t,W{t) + Un{t))dt — J^ g{t,W{t) + u{t)dt < 
(2M + l)e, and, since this holds for any e > 0, the result follows. □ 
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Note that Lemma [331 implies that pnk{x,y) and Pnk{x) are continuous, so that the esti- 
mates of Lemmas 13.11 and 13.21 will hold for all x,y & Q. 

We also need a stronger bound for sums of pnk terms than that given by the bounds for 
individual terms in Lemma [3.11 and the next two lemmas provide this. They are motivated 
by the idea that any solution of ([!]) should satisfy the approximate equation u{{k + 1)2"") ^ 
u{k2~^) + ijnk{u{k2~^)) which suggests that on a short time interval a solution can be 
approximated by an 'Euler scheme' Xk+i = Xk + <ynk{xk)- 

Lemma 3.5. Given even p > 2 we can find C > such that, for any choice of n,r G N 
with r < 2"/^, A; G {0, 1, ■ ■ ■ , 2" — r} and Xq G Q, if we define Xi, ■ ■ ■ ,Xr by the recurrence 
relation Xg+i = Xg + an,k+q{xq), then 

P (i2\Pn,k+q{xg^i,Xg) \ > 2-"|c5^|x,| + Ari/2|a;,||j < CX'P 
for any A > 0. 

Proof. We use Ci, ■ ■ ■ to denote constants which depend only on d and p. We write J-'j for 
T(k+j)2-"- Note first that Xq is jFg measurable and E(|cr„,/fc+g(xg )\P\Tq) < Ci2-''P/^\xq\P by 

dHD. Hence ¥.\an,k+q{xq)\P < Ci2~'''Pl^^Xq\'P . It follows that ^Xq+i]^ < {l + Cl^^2-''/yE\xq\P 
and so 

(16) E\xq\P < (1 + Cy^2-"/2)P|a;o|P < CsIxqI^ 

for 1 < q < r. 

Now let Yq = \pn,k+q{Xq^i,Xq)\, Zq = E(Yq\J='q) aud Xq = Yq - Zq. Then Xg Is Jg+i 

measurable and E(Xg|jFg) = so by Burkholder's inequality 



iK-i)r 



q=l 

< C,rP/^~-'2-^P/'J2^\xq - Xq_,\P = C5W2-i2-"P/2^E|(T„,fc+,_ 

r 

< C6rP/2"^2~"^'^E|x,_i|P 

9=1 

from which we deduce using f[TB|l that 

r 

(17) E| ^XJP < C7rP/22-"?'|a;o|P 

9=1 

Also let Vq = E{Zq\J^q^i) and Wq = Zq- Vq. Noting that Zq < C82-"/V„,g_i(xg_i) we get 
in a similar way that 

(18) E| J2 < C9rP/22-"P|xo|f 
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We also have 

(19) m < Cio2-^|X,_i| 

Now Yq = X, + Wq + Vq. By and (HID we have P(| Y!q=i{Xq + Wq)\> 2-''Xr^/^\xQ\) < 
Cn\P and the result then follows by □ 

Lemma 3.6. With probability 1 there exists C > such that for any n, r G N with r < 2"/^, 
any A; G {0, 1, ■ ■ ■ , 2" — r} and any yo, - ■ ■ ,yr & Q we have 

q=l \ 5=0 

where 7, = y^+i - yq - (Tn,k+q{yg)- 

Proof. Let Sn = 2~^"'^^. By Lemma [3.11 with probability 1 there exists C > such that, for 
any n,k > and any x,y & Q, we have 

(21) Pnk{x,y)<C2-^/'\x-y\+6n 

As before, let Qg = {x e M."^ : \\x\\ < 2'^}. Then, for integers s with < s < 2''/^, let 
be a set of not more than (2"(i^/^)^ points of Qs such that every x G is within distance 
q£ point of Qns and let fi„ = Uo<s<2-n/2f^ns- Let p = 8(4 + ci). Then by Lemma 1331 
there is Ci > such that the probability that 



r / r—l 

J2 \Pn,k+q{xq-uXq)\ > 2"" Ci ^ + X2^/'r'/' 

q=l \ q=0 



Xo\ 



for some n, r, as in the statement and some Xq G Qn, is bounded above by Ci ^'^^^ ^-P2"(3+d)2-pn/8 
which approaches as A — > oo. Hence with probability 1 there exists C > such that 



(22) Yl \pn,k+q{xq^i,x,)\ < C2-" Yl I'^il + 2"^'^'^' 

q=l \q=0 



r-1 

xo\ 



for all n, fc, r as above and Xq G f2ri. 

We now suppose, as we may with probability 1, that ( j2T]) and (!22l) hold (with the same 
C). We fix n, k,r,yQ ■ ■ ■ ?/r, 7o ■ ■ ■ 7r as in the statement of the lemma. Take the smallest 
s such that yo G Qs, noting that then 2^*^^ < ||/o| < (i^/^2^'^. Then we find Xq G fi„<( 
with |xo — yo\ < 2~'^~"' < 2-^~"||/o| and define Xi---Xr by the recurrence relation Xg^i = 
Xq + crn,k+qixq). Thcu by ([22D 



5^ |p„,fc+5(x,-i,x,)| < C2-" J2 + 2"^' 

g=l \q=0 



r-1 

I I . ^n/4l I 
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Using ([21]) we have = \xg + an,k+qi^q) \ ^ (1 + + 5„ so \xg\ < Ci(|xo| +r5„) 

and 

r 

(23) 5^|p„,fc+,(x,_i,x,)| < C22-=^"/^(|xo| +2"/^) 

9=1 

Now let Ug = Xg- yg. Then - Ug\ < \pn,k+q{xg, yg) \ + |7g| so 

< |m,|(1 + C2'"/^) + |7,|+5„ 
and since |uo| < 2^~'*|?/o| we deduce that \ug\ < C3(2^"|?/o| + r5n + Xlg=o l7gl) ^^"^ so 

(r-l 
2^"bo|+r5„ + 5^|7,| 

and we have the same bound for |p„_fc+g(xg_i, yg_i)|. Now 

Pn,k+g{yg-l, Vq) = Pn,k+q{Xg^i, Xg) + Pn,k+q{yq-l, Xg^i) + Pn,k+q{Xg, yg) 

and then using (!23|l . (l2ll) and the fact that |xo — l/o| ^ 2^^"|?/o| we deduce that 

r / r—1 ^ 

\PnMq{yq-u yq\ < c, 2-=^"/^|yo| + 2-"/^ ^ |7, 1 + 2-"/25„ 

g=l V 9=0 

from which the result follows. □ 

We now proceed to complete the proof of the theorem. From now on we take g = f in 
the definition of ank and pnk- We consider a Brownian path W satisfying the conclusions of 
Lemmas 13.11 13.21 13.61 and 13.41 for some C > 0. We shall show that for such a Brownian path 
the only solution m of (jlj) in $ is m = 0. This will follow from the following: 

Lemma 3.7. Suppose W satisfies the conclusions of Lemmas \3.1[ \3.S\. \3.6\ and \3.4\ for some 
C > 0. Then there are positive constants K and rriQ such that, for all integers m > mo, if u 
is a solution o/ m $ and for some j G {0, 1, ■ ■ ■ , 2'" — 1 and some /3 with 2~2^™'''' < (3 < 
2-2''"/' we have |m(j2-'")| < p, then 

\u{{j + 1)2'"^)\ < f3{l + K2~"^ log(l//5)} 

Proof. We use Ci,C2, ■ ■ ■ for positive constants which depend only on the constant C and 
the dimension d. Fix m, j and /3 as in the statement, and suppose \u{j2~"^)\ < [3. Let N be 
the integer part of 41og2(l//3). Suppose m G $ satisfies ([1]), and let m„ be the step function 
which takes the constant value u{k2~^) on the interval /„fc, for = 0, 1, ■ ■ ■ , 2" — 1. 
Let a be the smallest nonnegative number such that 

(j+l)2"-'"-l 

(25) \u{{k + 1)2-") - m(A;2-") | < a2-'''\n^'^T''^ + N) 

for all n with m < n < N . 
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For n > m let 

(j+l)2"-'"-l 

(26) V'n= Yl l«(^2-")| 
Then by ([25]) 

for n > m,and since ipm = /? it follows that 

n 

(27) < 2"-"^/? + ^ «2"~'~™(/^/22V2 + A^) < Ci2"""^(/? + a2-'^N) 

l=m+l 

for all n with m < n < N, where we have used the fact that 

^l/22m/2 jg 

bounded by const. A^. 
Now fix n > m. Then for k = j2"~"^, ■ ■ ■ + 1)2"~™ — 1 we have, using ( |T5l) 



«((fc + l)2-")-n(A:2-")= / {/(iy(t) + ^t)) - /(iy(t))}rft 

' {f{W{t)+Ur,it))-f{Wmdt + Y, / {/(H^(t)+«mW)-/Wt)+«^W)}^i 



l=n 

which we can write as 
(28) 



oo (fc+l)2'-"-l 



u{{k + 1)2-") - u{k2-^) = ank{u{k2-'')) + ^ pi+^^2r+iH2-'-\2r + 1)), m(2-V)) 

from which we deduce 

(i+l)2"-'"-l 0-+l)2"-'"-l oo 

(29) 5^ |^i((A; + l)2-")-M(A;2-")| < |a„fc(M(A;2-"))| + ^ fi^ 

k=j2"-"^ k=j2"-"^ l=n 

where fi, = ^15'/-^^' |p;+i,2.+i(M(2-'-H2r + 1)), m(2-V))|. 

We now proceed to estimate the two sums on the right of (129!) . starting with the easier ank 
term. Using LemmaOand the fact that < 2"", we have |a„fc(x)| < C2n^^'^2-''/'^{2~^ + \x\) 
and so 

(i+l)2"-'"-l (j+l)2"-"'-l 

E l^"^(^(^2-"))| <C, ^'/'2-/2(2-^ + |^(A:2-")|) 

< C3n^/22n/2-m^^ ^ 2-^iVa + 2-^) 

using (1271). 

Next we bound X]^'> which we do in two stages. We first obtain a relatively crude 
bound by applying (fTTj) to each term, and then obtain an improved by applying the crude 
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bound together with Lemma (13. 6p . To start with the crude bound, from ( |TT1) we have 
\Pnkix,y)\ < C32""/^A^^/^(2~^ + \x — y\) and using this together with fl25l) gives 

(31) ill < C42"'/2Ari/2{2-^2'-'" + a2~"'{l'/^2^/^ + N)} 

and so 

N 



(32) XI - ^5 ^p^l/22-m-N/2 ^ c,2-"'N^) 

l=m 

For I > N we use |-u(t) - ^(t')! < |t - t'\ and ([11]) to obtain 

oo oo 

(33) Yl C'62'-'"/'/'2-3'/2 < C7iVV22— 

Z=Ar+l i=Ar+l 

and combining this with (1321) we obtain 

oo 

(34) Y.^^- C8(A^'/'2-'"-^/2 ^ a2-"'N^) 

l=m 

The second stage is to improve the estimate ( IMI) by applying Lemma 13.61 to obtain a 
better estimate for Qn for larger n; we use ([31]) to bound the 7 term in Lemma [3.61 

Let N^/^ <n<N. We define 7^^ = m((A; + 1)2"") - u{k2-'') - a„fc(u(A;2-")), noting that 
fl28]) imphes that 

(i+l)2"-'"-l 00 

(35) X |7„fc| <Y^i< C8(iVi/22— + a2-™iV2) 

fc=j2"-™ l=n 

Also we define 

(i+i)2"-™-i 
An= J] |p„,fe+i(2-"A;,2-"(A; + l)) 

SO that ^In < A„_|_i. Let r = [2"/^J. In order to apply Lemma [3.61 to estimate A„, we will 
split the sum into r-sized pieces. First we find i G {0, 1, ■ ■ ■ , r — 1} such that, writing s = 
[r^^ {2"'^'"^ — i) \ , we have Ylt=o\'^W^~'^ ~^ ~^^^)'^^^)\ — '^~^'^n- Now we fix for the moment 
t e {0, 1, ■ ■ ■ , s} and apply Lemma 13.61 with yq = u{{k + g)2~") where k = j2^~"^ + i + tr. 
We obtain 

r / r—1 

J2 \PnM,iy,~i,y,)\ < C9 2-3"/V(A;2-")| + 2-"/^X |7.,.+,| + 2-'"^' 

9=1 \ q=0 
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Summing over t then gives 

(i+l)2"-'"-l s 

J2 |Pn,A:+i(2-"A;, 2-"(fc + 1))| <C,2-'^/' Yl + + ^02"")| 

((j+l)2"-'"-l 
+2--/^ 5^ |7n,.|+2"-^"^^ 

Also 

i2"-'"+i-l / j2"-™+j-l 

fc=j2"-'" \ fc=j2"-™ 

From the last two inequalities, using (1271) . (!35|) and |'u(j2^"*)| < /?, we find that 

An < Cio{2-'"(/3 + a2-'"Ar) + 2-^-/4 (AT V22~^/2 + ^tAT^) + 2"-2"''^} 

bmce n > N^/^ the first term dominates so A„ < Cii2-"^{(3 + a2-'^N), and the same bound 
holds for Vtn < An+i- We deduce that 

m/(i<i<N 

Using the original bound fl?I]) for I < N^/^ we have 

Ci3ArV2(2-^+^''V2-m + a2-'^{N^I^ + 2-^/''N)} 

Combining these two estimates with fl33|) we get our improved bound. 



l=m 

To conclude the proof we use this bound along with fl30l) in fl29l) and obtain 

(i+i)2"-'"-i 

^ \u{{k + 1)2-") - u{k2-")\ <Ci5(ni/22"/2-™ + iV2-"^) 

X {/5 + a{N2-'^ + N-^'^ + 2-'^/2Ari/2)} 
for all n with m < n < N . Comparing this with (1251) we see by the minimality of a that 

a < Ci5{/3 + a{N2-"' + N~^'^ + 2-'"'^N^'^)] 

Then if m is large enough to ensure Ci^{N2-"^ + N^^'^ + 2-™/2jyi/2) ^ it follows that 
a < 2Ci5/?. Then applying ([25]) with n = m gives + 1)2""^) | < /? + 2Ci5/3(mV22™/2 + 
A^)2-™ < /3(1 + Ci6A^2-'^) from which the required result follows. □ 
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To complete the proof of Theorem 11.11 using the notation of Lemma 13.71 let m > itlq 
and f3o = 2"^^"^*, and define Pj for j = 1,2, ■■■ ,2™ by the recurrence relation /3j_^_i = 
f3j{l + K2-"^log{l/(3j)). Writing 7^- = log(l//?j) we then have 

7,+i = 7. - log(l + K2-7,) > 7,(1 - K2~"^) 

so the sequence (7,) is decreasing and 

Ij > 7o(l - K2~''^y > -foe-^-^ = 23W4g-x-i > 22^/3 

for all J = 1,2, ■ ■ ■ ,2™, provided m is large enough. Then for each j, f3j is in the range 
specified in Lemma 13771 and it follows from that lemma by induction on j that |M(j2~™) | < fij 
for each j. Hence \u{j2~'^)\ < 2^^^™ ^ for each j. This holds for all large enough m, and 
hence u vanishes at all dyadic points in [0,1], and, as u is continuous, m = on [0,1]. This 
completes the proof of the theorem. 



4. An Application 

We give an application of Theorem II. II to convergence of Euler approximations to ([1]) with 
variable step size. 

In this section we assume / is continuous and consider ([1]) on a bounded interval [0,T]. 
Given a partition V = {0 = to < ti < ■ ■ ■ < = T} of [0,T] we consider the Euler 
approximation to ([1]) given by: 



for n = 0, ■ ■ ■ , — 1, with xq = 0. For such a partition V we let S(V) = max^^^(t„ — tn-i) 



Then we have the following: 

Corollary 4.1. For almost every Brownian path W , for any sequence 

n = {ti'\---,t^Pj 

of partitions with S{Vk) 0, we have 



max 

n=l 



as k oo, where x{t) is the unique solution of (QP and {x^f^} is the Euler approximation 
using the partition Vk- 

Proof. Suppose is a path for which the conclusion of Theorem 1 1 . 1 1 holds . and suppose there 
is a sequence of partitions with 6{Vk) —>■ such that max^^^^ — x{tn ^)\ > 6 > 0. Then if 
we let Un'^ = Xn^ — W{ti!''^) we have \ul^li—Un^\ < ||/||oo(^i+i — ^ri^^) so by Ascoli-Arzela, after 
passing to a subsequence we have a continuous u on [0, T] such that max^^;^ \un^ — u{tn^)\ — > 
0. Then writing y(t) = u(t) + W(t) we see that y ^ x and, using the continuity of /, that y 
satisfies ([1]), contradicting the conclusion of the theorem. Corollary 14. II is proved. □ 
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The point of Corollary 14.11 is that the partitions can be chosen arbitrarily, no 'non- 
anticipating' condition is required. For general SDE's with non-additive noise and suffi- 
ciently smooth coefficients Euler approximations will converge to the solution provided the 
partition points t„ are stopping times, but this condition is rather restrictive for numerical 
practice, and an example is given in section 4.1 of [1] of a natural variable step-size Euler 
scheme for a simple SDE which converges to the wrong limit. [1] also contains related results 
and discussion. 
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